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Abstract. Explicit formulae are given for the Airy and Bessel bispectral involu- 
tions, in terms of Calogero-Moser pairs. Hamiltonian structure of the motion of 
the poles of the operators is discussed. 



1. Introduction 

This paper follows upon the study of the Airy bispectral involution made in j[KR . 



There we gave an analogue, for arbitrary rank, of the rank-one bispectral involution 
developed by Wilson ||W1|| . Recently, ||W2j"| , Wilson has established a relationship 



between the rank-one bispectral involution and the complex analogue of the Calogero- 
Moser phase space. This relationship leads to explicit formulae for the Baker function 
and the corresponding involution, which make many important features manifest. As 
shown below, similar results hold for bispectral algebras obtained from generalized 
Airy and Bessel operators. 

Given a positive integer n, define C n to be the quotient, under conjugation by 
Gl(n,C), of the space of pairs (Calogero-Moser pairs) of n x n complex matrices, 
(P,Q), such that 

(1.1) rank([P, Q] - I) = 1 . 



This is the complex analogue of the definition in |KK5 |, in which P and Q are taken 



to be hermitian, and / is replaced by il. Define Gr ad (see ||W 1| , |W2[ ) to be the 
subspace of the Sato grassmannian ([ pjKM| , |5a| , |S W| |) corresponding to Krichever 



data whose spectral curve is rational and unicursal (no nodes). The Baker function 
for such data is always of the form e xz p(x, z), where p(x, z) is rational and separable. 
By separable we mean its denominator is a product q(z)r(x). Wilson has proved 

Theorem 1.1 ( ||W2|| , Thm. 5.1). There is a one-to-one correspondence 

(1.2) U™ =0 C n ^Gr ad , 
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such that a point W G Gr ad corresponds to a point (P, Q) if and only if its Baker 
function, ipw is given by 



(1.3) ifj w (x, z) = e xz Det(I - (zl - Q)~\xl - P)" 1 ) . 

Remark: The "generic" rational unicursal curve has simple cusps. These corre- 
spond to Q with distinct eigenvalues. In the hermitian case this is automatic. An 
important aspect of Wilson's theorem is the observation that in the complex case, 
nonsemisimple Q correspond to nonsimple cusps, or in physical terms, collision of the 
Calogero-Moser particles, which are now moving in C rather than BL 

The spectral algebra IZw is the algebra of differential operators L(x, d) such that 
Lijj w = f(z)i/jw for some function f(z). Wilson's result proves that the spectral 
algebra of any point W G Gr ad is bispectral in the sense of |pG| . That is, the 
spectral algebra of ipw(z,x) is also nontrivial. 

It proves much more, for it says that tpw(z,x) is also a Baker function, namely 
the one corresponding to (Q T ,P T ). The involution (P,Q) i— > (Q T , P T ) is clearly 
antisymplectic with respect to the symplectic form 

(1.4) lo = tr(dP A dQ) . 

This symplectic structure is an important example of "unreduction" [|KKS|| . Namely, 
the Calogero-Moser hierarchy ||AMM|| is a completely integrable hamiltonian system 



defined on the quotient space C n . The hamiltonians are rather complicated in the 
reduced variables, but on the level of matrices (P, Q) they are given simply by the 
hamiltonians are h n = tr(P n ). Moreover, the involution ((P, Q)) — > ((Q T ,P T )) 
is the linearizing map for the Calogero-Moser particle system ||AMM| , [KKSf . Thus, 
Wilson's result gives the best proof that this linearizing map and the bispectral 
involution are one and the same (cf. |[Ka|| ). As we shall see, the corresponding 
involutions, in terms of an auxilliary monic polynomial p(t), are given in the Airy 
and Bessel cases respectively by 

(1.5) (P, Q) ^ (P T , p(P T ) - Q T ) , 

(1.6) (P, Q) » {{QPp{QP)- l Q) T , (Q- 1 p(QP)) T ) . 



Acknowledgments: The author wishes to acknowledge valuable discussions with J. 
Harnad, A. Kasman, R. Varley and G. Wilson. 



2. The Airy case 



The rank of a commutative algebra 1Z of ordinary differential operators is defined 
to be the greatest common divisor of the orders of its elements. The true rank of 1Z is 



BISPECTRAL INVOLUTIONS, February 9, 2008 3 

the rank of its centralizer |[LP| , [P W|] . For instance, fix a monic rth order polynomial 
p(t), and define 

(2.1) L Ai = p(d) - x (Airy) 



Then C[Z/Ai] has true rank r [|KR|| . Moreover, this algebra is bispectral. Indeed, for 
any / G ker(L A {), define 

(2.2) f Ai (x,z) = f(x + z) . 
Then 

(2.3) L Ai (f Ai ) = zf Ai , L Ai (z,d z )(f Ai ) = xf Ai . 

To generalize Wilson's formula to this case, introduce a Baker functional, ^m,Wj 
defined on kerL Ai , of the form 

(2.4) *Ai(/) = E**fo • 

8=0 

To specify its properties, it is useful to introduce the dual description of Gr ad (cf. 

As Wilson shows, every point of W G Q r ad arises from a homogeneous, 



finite dimensional space of finitely supported distributions in the complex plane. That 
is, there should be complex numbers Ai, . . . , A„ and polynomials over C, l\, . . . , £ n , 
such that 

(2.5) W = ^-r { p(z) e C[z] \ Cl (p) = 0} , 
where 

(2.6) q{z) = U i {z-\ i ) 
and 

(2.7) a = 5 Al o t{d z ) . 

Now define ^ A \ t w by the following properties. Let C denote the span of c i; c n . 
Property la: The functions q(z)ki( polynomial in z. 

Property 2a: For all / G Ker(L), qc(z)^> 'M,w(f) is annihilated by all c G C. 
Property 3a: lim (ho, fc r _i) = (1, 0, 0). 
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Fixing the operator L A {, set 



k 



(2.8) 



k 



r—l , 



It is important to remark that the map C \— > W is not one-to-one. The equivalence 
relation on C induced by this map is the one generated by 

(2.9) C + C8 X ~ Co (z-X) . 

In particular, the properties of ^Ai,w depend only on W and not on the representative 
C. 

Choose a space of conditions C defining W, and set ki(x,z) = ■^j^J2jh,j( x ) z: ' 
Consider the differential operator 

(2.10) K Ai:C = Y,hj( x )9 l L 



It TO 
Ai 



'■J 



Defining 

(2.H) fAi(x,z) = (/Ai^/Ai),...,^ 1 ^)) , 

one has 

(2.12) K AUC (f Ai )=q c (z)f M -k w . 

The asympotics of kw imply that i^Ai,c is a monic operator of order rn. Property 
2a implies that K\^c annihilates c(/ai) for all the distributions c G C and all / G 
Ker(L A {). Thus K A i,c is unique operator with the two properties just stated. It then 
follows (cf. ||KR|1 ) that for any polynomial p G Rw, the pseudo differential operator 

c is a differential operator, and 

(2.13) M p (K Ai ,c(f M )) = p(z)K AiyC (f Ai) 

for all / G Ker(L A i). Note that M p could also have been obtained by conjugating 
p(L Ai ) by the monic 0th order pseudodifferential operator K A ^ c <i{.L Ai )~ l . The latter 
operator is independent of the space C representing W, and is the analogue of the 
Sato operator in this theory. 

Thus one has a rank r commutative algebra of differential operators 

(2.14) K M , W = {M p \ P e R w } , 

with an r-dimensional space of eigenf unctions ty A i t w{f)> f £ Ker(L). To accomodate 
the usual normalization of spectral algebras, namely that the subprincipal symbol 
should vanish, take p of the form p(t) = t r + 0(t r ~ 2 ). 
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The task is to obtain a formula for kw in terms of the matrices P and Q corre- 
sponding to the point W. The following lemma expresses property 2a in terms of 
covariant differentiation of k w . Set 



(2.15) 



B Ai (x, z) 



x + z 

a r _ 2 
1 



where the a^'s are the coefficients of p, and set 

d 



(2.16) 



V A i = ^ + B Ai (x, z) 



9o(x,z) 

Lemma 2.1. Let c be a distribution of the form c — 8\ op(d z ). Let g = 

_g r -i(x,z) 

— * 

be a vector of polynomials in z with coefficients in C(x). Then c{fA% • g) — for all 
f G Ker(L) if and only if g is annihilated by 5\ o p(V Ai)- 

Proof. Given / e Ker(L), 

(2.17) d z (f Ai ) = f A i-B Ai (x,z) . 
Thus, for all j, 

(2.18) S x o di(f Ai ■ g) = f Ai (x, A) • (d z + B Ai (x, . 

This proves the lemma, since there is no differential equation of order less than r 
satisfied by f(x + A) for all / G Ker(L). □ 

Now consider the involutions on each C n defined by 

n. ^n. 



(2.19) 



((P,Q))^((P,Q)) , where 
P = P T - Q = P (P T )-Q T . 
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Theorem 2.2. Let W G Gr ad correspond to a point ((P, Q)) G C n . Let 



(2.20) 



[P, Q]=I- w lW \ 



where w\ and W2 are column vectors. For j = 0, . . . , r — 1, let 
(2.21) 



r-2 



Pj (t)=f- i -j- y: att 

i=j+l 

— * 

Then the components of k w are 

(2.22) kj(x, z) = 5 ,j - w T 2 {zI - Q)- 1 p j (P)(x - Q T )~ 1 w 1 . 

Proof. It suffices to consider the generic case, in which W is defined by conditions 
Ci — S\ { o (d — cxi), i — 1, • • • , n, for distinct A's. Set 



(2.23) 



7» = «< - S 



/ • A, — Ao 



Then corresponds under Wilson's theorem to the Calogero-Moser pair 

"Ai 

(2.24) Q 







At 



(2.25) 



7i 



A2— Ai 



1 

An— Al 



1 



Ai— A2 



Al -\n 

A2— A n 



In 



Let (ej)j = i t ,_, tr be the standard basis for C r . Then 



(2.26) 



k w = ei + Yl 



VAX 



for some vectors Vi(x) G C r . Applying 5\ % o (d z + E>Ai(x,\i) — ati) to q(z)kw, one 
obtains the set of equations 

= B Ai (x, \)vi(x) H(\ - X e ) + 

n 

n(A, - A,) ei + v m {x) e n ( x i - x s) 

(2.27) - ai v, t (x)Y[(\ t -\ E ), 
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i — 1, . . . , n. Dividing by Yie.^i{\ ~~ ^t) anc ^ usm g (2-23), 

v e (x) 



(2.28) 



ex = -B Ai (x, Xi)vi(x) + 'jiVi(x) -J2 \ 

e^i A i ~ ^ 



Let ( 

Ui)i=i,...,n be the standard basis for C n , and let 
(2.29) v(x) = J2ui®Vi(x) eC n ®C r . 

Let w = J2 u i- Write B Ai (x, Aj) = B Ai (x, 0) + AjA, where 



/0 




1\ 




(2.30) A = 

\0 ••• 0/ 
Then (2.28) is encoded as a single equation 

(2.31) w <g> ei = (-J <g> B(x) + P ® I - Q® A)v(x). 
Altogether, (2.26) becomes 

(2.32) k w = e 1 - (w T <g> I) o ((zl - Q)' 1 <g> I) o A~ 1 {w <g> a), 
where 

(2.33) A = I <g> 5 Ai (z, 0)-P®/ + Q®Ae £nd(C™ <g> C r ). 
Thinking of C n ® C r as rn-tuples in blocks of length n, 

"11 



(2.34) 



w 05 ei 



>n 



and 



(2.35) 



-P ••• + 
J -P aj 

o / ••. ; 



o 





a r -\l 
I -P 
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One checks quite easily that 



(2.36) 



Then 

{w T <g> I) o ((zl - Q)" 1 ®I)o A~ 1 (w ® d) 
,t 



" Po(P) ' 




'xi + q- P (py 









Pr-l(P)_ 








w 







w 



zI-Q)' 1 
(zl - Q)- 1 



Po(P) 



Pr-l(P). 



(x/ + Q-p(P))-V 



(2.37) 



w T (zI - Q)- 1 Po (P)(xI - Q T )~ l w 

_W T (ZI - g)- 1 pr-l( J P)(x/ - Q T )~ 1 W_ 

This proves the theorem, since 
(2.38) [P, Q] = I- ww T . 

□ 



Corollary 2.3. For any generalized Airy operator La%, and for all W G Gr ad , 

(2.39) k w (z,x) =k MW) (x,z) . 

In particular, the algebra TZai,w, (2-14) ; is bispectral, with an r- dimensional space of 
joint eigenf unctions \1/ Ai,w(f), f £ Ker(LAi). 

Proof. Formula (2.39) follows immediately from (2.22), for if [P, Q] = I — wiw^, then 

(2.40) [P, Q] = [P T , p{P T ) - Q T ] = [P, Q] T = I- w 2 wj . 
The rest of the corollary is immediate. □ 
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3. The Bessel Case 

The Bessel case works in much the same way. Consider again a polynomial p(t), 
now normalized so that a r _i = Q • Set 



(3.1) 
where 

(3.2) 



L Be = x r p(D) (Bessel) , 



d = — , D — xd . 

ax 



Consider Ker(L Bc — 1), which should now be thought of as a sheaf rather than a 
space. For / e Ker(L Bc — 1), define 

(3.3) /Befa,*) = /(xz) • 
Then 

(3.4) L Bc (f Bc ) = z r f Bc , L Bc (z,d z )(f Bc ) = x r f Be . 

Assume now that the matrix Q is invertible. Define, as the analogue of (2.11), 

(3.5) 
Then 
(3.6) 



/^(x^) = (/Be, J D(/Be),..., J D r - 1 (/Be)) 

D z (f Bc ) = f Bc B Bc (x r ,z r ) , 



where 



(3.7) 



B Be (x, u) 



... «o + xu 

1 ' • . Oi 
'•• : a 2 

a r _ 2 

... 1 a r _! 



Accordingly, one expects a Baker functional of the form 

(3.8) *Be,w(f) = /Be " ^) , / G ^er(L Bc - 1) • 
To state the properties of ^Be,^, introduce the functions 

(3.9) fi(x,z) = (x r ,z r ) ; u(z) = z r . 

Denote by v* the action of v on the space of finitely supported distributions in C*. 
Let 



(3.10) 



V Be = D z + -B Be (x, z) . 
r 
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Given a distribution c = 5\ op(D), define 

(3.11) c Vbc = 6\ op(V Be ) , 
acting on vector valued functions of z. 

Lemma 3.1. Let c be a distribution of the form c = 5\ op(D). Let g = 

g r -i(x,z) 

be a vector of polynomials in z with coefficients in C(x). Then cyfse 4 1 1 *^)) = f or 
all f G KeriLiBe — 1) if and only if v* : (c)v Be (g) — 0. 

Proof. By virtue of the identity 

(3.12) Dou* = u*orD, 
one has 

(3.13) u*{c) = 8 X rop(rD) . 
By (3.6), 

co f Bc - n*{g) = 5 X o f Bc ■ p(D + B Be (x r , z r )) o ^*(g) 

= fBe(x, A) ■ S x o /i* op(rD + B Bc (x,z))(g) 
= fBe(x, A) • v* x o 5 X r op(r\7 Be )(g) 

(3.14) =f Be (x,\)-v*ov*(c) WB M , 



where u x is v acting in the x-variable. The lemma now follows as in lemma 2.1. □ 



In light of the preceding lemma, it makes sense to impose the following properties 

On ^Be,W- 

Property lb: The functions qc(z)ki( 

X) 21X6 polynomial in z. 

Property 2b: Let C be any space of distributions such that u*(C) = C. Then for 
all / G Ker(L Be — 1), qc( z )^Bc,w(f) is annihilated by all c G C 

Property 3b: lim k w = e\. 

As in the Airy case, one reconstructs a differential operator K Be C , but now 

(3.15) #Be,c(/Be) = qcffifee ' k W (x T , Z r ) . 
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Then for any polynomial p G R\y, the pseudo differential operator M p = K-Q e) cp{L-Q e )K B l c 
is a differential operator, and 



(3.16) 



M p (K Be)C (f Be ))=p(z r )K Be;C (fB e 



(3.21) 



= 5 Xl o (d z + — B Be (x, Xt) - ai)(qc{z)k w (x,z)) . 



for all / G Ker(L Be — 1). Define 72-Be,w to be the algebra of the all the M p 's. 

Everything now proceeds as before. Assume that the matrix Q is invertible. We 
have n distributions q = 5\ t o (d z — aij). Note that 5\ i o (d z — cuj) = j-Sx, ° (D z — AjO;,). 
Thus, according to lemma |3.1| , property 2b imposes the n conditions 



(3.17) 
Setting 

(3.18) 

one now finds 
(3.19) 

This time, 

(3.20) 

where 



" Vi{x) 



ei 



•-^-5 Bc (x, Ai)tfi(x) + 7i?Ti(x) - ^ r^—\ ■ 

T"Ai ^ Ai — A£ 



1 1 X 

—B Be (x,u) = —Ax + ~A 2 , 



ru 



ru 



Ax 



a 



Ola 



r-l 



(3.22) 

Thus 
(3.23) 



A, 



o ••■ o i 
■ ■ o 

o ••• o o 



x 1 
w <g> e 1 = -(-I ® A 2 - P <g> J + -Q" 1 ® Ax)u(x) 
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Then 

(3.24) k w = e x - (w T ® I) o ((zl - Q)' 1 ® I) o A~ l (w <g> ei), 

where A is now given in block matrix form by 

-P ••■ + 



(3.25) 



.4 



r ^ 










fQ- 1 

a.r-2 D—l 



-p 

IQ- 1 -P + V^ 1 



One obtains the following result. 

Theorem 3.2. Let W E Gr ad correspond to a point ((P,Q)) E C n . Let 
(3.26) [P, Q]=I- wiv% 

where w± and u> 2 are column vectors. Writing the rth order Bessel operator L# e in 
the form L Be = x~ r p{D), let 



(3.27) 



Pj (t)=t' 



i=j+i 



for j = 0, . . . , r — 1. T/ien £/ie components of kw are 

(3.28) fcj-fo z) = 5o d - ru%{zl - Q)~ l P] (rQP)(x - Q T )~V ; 
where 

(3.29) Q = {Q- l p{rQP)) T . 
Proof. One checks now that with A given by (3.25), 



(3.30) 



.4 



rp (rQP) 
rpr^rQP) 



xl - Q- x p(rQP) 




The result then follows as in theorem 12.21. □ 



Theorem |3.2| suggests the definition 

P = Q- 1 P T Q T 
(3.31) = (QPp(QP)- l Q) T 
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Then 



Q = (Q-V(rQP)) T 

= (g-V(rP T Q T )) T 

(3.32) = p(rQP)p(rQP)- 1 Q = Q 



and 



(3.33) = Q~ l QP{Q~ l ) T Q T = P ■ 
Moreover, 

[P,Q] = [Q- l P T Q T ,Q] 

= Q~ 1 P T Q T Q - P T Q T 

= Q T {p{rQPY l ) T {QP) T p{rQP) T {Q- l ) T - P T Q T 

(3.34) = Q T P T - P T Q T = [P, Qf . 

So again one has an involution (densely defined) on C n , 

c --Uc 

(P, Q) i > {P, Q) , where 

(3.35) P = Q- 1 P T Q T ; Q = (Q" 1 p(rQP)) T . 



Corollary 3.3. For any generalized Bessel operator Lb 6 , 
(3.36) k w (z,x) = kp Be{w) (x,z) . 

In particular, the algebralZBew is bispectral, with a rank r joint eigensheaf^Bew(f), 
fe Ker(L Be -l). 

4. Dynamics 

It is well-known that the Calogero-Moser particle system is a completely integrable 
hamiltonian system on the symplectic manifold C n . The symplectic form is 

(4.1) u = ti(dPalQ) , 

and the hamiltonians are h n = tr(P n ) (cf. |[KKS|| ). 

It is pleasing then that the Airy and Bessel involutions are antisymplectic on each 
C n . In the Airy case, 

(4.2) P = P T ; Q = p (P T )-Q T , 
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this follows from the fact that 

(4.3) tr(dPdP n ) = ^(dPFdPPi) = . 

i+j=n-l 

(The basic trace identity for form- valued matrices is 

(4.4) tr(XY) = (-l) dcs W dcg ( Y Hr(FX) .) 
The Bessel case is slightly more involved. 

Proposition 4.1. Let a be a polynomial over C. Let 

(4.5) Q = (Q- 1 a(QP)) T ; P = Q- 1 P T Q T . 
Then 

(4.6) tr(dPdQ) = -tr(dPdQ) . 
Proof. Let Q = Q-^QP), P = QPa{QP)~ 1 Q = P T . Then 

(4.7) tr(dPdQ) = tr(dPdQ) . 
Set R = QP. Then 

(4.8) tr(dPdQ) = d tr(Q _1 i2dQ) , 
while 

tr(dPdQ) = dtr(RxT- 1 Qd(Q- 1 a)) 

= d tY{Ra- l da - Ra^dQQ^a) 

(4.9) = d t^Ra^da - Q' l RdQ) . 

So it must be proved that 

(4.10) dtriRa^da) = . 
If a = <J\(j2 and o 2 commutes with R, then 

d ti(Ra~ l da) = dtT(Ra 2 1 a^ 1 (da 1 a 2 + oidc^)) 

(4.11) = d t^Ra^da^ + d t^ifcr^dos) . 
Also, 

(4.12) d tr(a -1 d<T) = -tr^'Maa -1 ^) = 0. 

This last identity implies that one can replace R by R— (const) I in (4.10), and reduce 
to the case that a = Rai(R), o\ polynomial. 
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Then by (4.11) 

d ti(Ra' l da) = d tr{dR + Ra^ l da{) 

(4.13) =dtr{R<r± 1 da 1 ) . 

Since the result is obvious when a is a constant, the proposition follows by induction 
on the degree of a. □ 

Now introduce time dependence into the Baker functionals of the previous sections 
in a manner generalizing the standard procedure in the rank-one case [ |SW|| . Fix a 
positive integer m. The time-dependent Baker function in the rank one case is the 
function 

(4.14) e xz+tzm p(x,z, t) 

satisfying properties 1,2 and 3 of section |2| with a fixed space of conditions C. On 
the other hand, one may introduce time-dependence into the conditions by defining 
Ct = C o e tzm . Then the function e xz p(x, z,t) satisfies properties 1,2 and 3 for the 
variable conditions C t . 

If c = 5\ o (d — a), then c o e* 2 ™ = 5\ ° (d + tm\ m ~ 1 — a). In other words, the flow 
Ct is seen on the level Calogero-Moser pairs as 

(4.15) Q t = Q ; P t = P - tmQ m ~ l . 
This is the flow of the completely integrable hamiltonion 

(4.16) hm = tr(Q m ) . 

This hamiltonian is the Calogero-Moser hamiltonian with the roles of P and Q re- 
versed. Finally, the Baker function gives rise in a standard way to a solution of the 
KP-hierarchy, with poles in x the same as those of the Baker function. Thus, Wil- 
son's formula makes it immediately clear that the poles in x of such a KP-solution 
move as a Calogero-Moser particle system (cf. ||Kr| , f5h] , [Ka|| ) . 

To carry this over to the Airy and Bessel cases, define k Wjt to be the vector (2.8), 
for the variable space of conditions C t . We are led to solutions of a subhierarchy of 
the KP-heirarchy, in the following way. Let i^o.Ai be the monic 0* fe -order pseudodif- 
ferential operator such that 

(4.17) K 0Ai d r K^ M = L Ai . 

Let K t be the monic 0^-order pseudodifferential operator such that 

(4.18) fM-kw,t = Kt(fAi) ■ 
Then the argument in [[5W11 shows that 

(4.19) d t (k)R- 1 + {KLZK- 1 )- = . 
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Now let 

(4.20) K t = K t K 0>M . 
Then 

(4.21) d t (K)K- 1 + (Kd rm K~ l )_ = . 
It now follows as in [ R W|| that the operator 

(4.22) M t = KtdK- 1 
satisfies the rm th term of the KP-hierarchy, 

(4.23) d t (M) = [M r + m ,M] . 

From formula (2.22), M has coefficients in C(t)[x, ^r^y]? where 

(4.24) Tt (x) = Det(x - Q t ) , 



Q t = (p(P t ) - Q 



(4.25) =(p(Po-tmQ^ L )-QoV ■ 

Thus we have constructed a solution M t , of the rm th term of the KP-hierarchy, whose 
poles in x move according to the completely integrable hamiltonian tr(Q m ) on the 
Calogero-Moser phase space. 

In the Bessel case, exactly the same analysis holds, except that M has coefficients 
in C(t)[x, |, Tt Lr\ ]- Thus M t has a fixed pole at x — 0, with the motion of the 

remaining poles being governed by the hamiltonian ti{Q m ). 

With several particles it becomes quite cumbersome to write out these hamiltonians 
explicitly. The lowest rank cases are 

(4.26) L Ai = d 2 - x ; p(t) = t 2 , 

(4.27) L Be = d 2 - x~ 2 ; p(t) = t 2 - t - 1 . 

Note that the first hamiltonian, h\ = tr(Q), is already non-linear. In the rank-two 
Bessel case with one particle, for instance, h\(\, 7) = A _1 p(A7). This gives the 
equations of motion 

(4.28) A = 87A - 2 

(4.29) 7 = -4 7 2 - A~ 2 . 

These equations are solved by applying the Bessel involution and changing t to —t, 
i.e. by setting 

(4.30) A = ci ; 7 = c 2 + t . 
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X(t) = 4cit 2 + (8cic 2 -2)t 2c 2 + 4ciCn . 

Cl 



After some calculation, 
(4.31) 

With two particles and a second order Airy operator, the first hamiltonian is 
(4-32) 7l 2 + 72 2 - Ai - A 2 - - 2 . 

{ — Ai + A 2 J 

With two particles and a second order Bessel operator, the first hamiltonian is 



(4.33) 

[AMM] 

[BHY1] 

[BHY2] 

[BHY3] 

[DG] 

[DJKM] 

[Ka] 

[KKS] 

[Kr] 

[KR] 
[LP] 

[PW] 

[Sa] 

[Sh] 



Ai + A 2 . 2 , * 2 2 ( A i 7i - A 2 72) 

- 7i + M 7i - 72 + A 2 72 + - 



Ax A 2 



-Ai + A 2 
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